Electron delocalization and multifractal scaling in electrified random chains 
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Electron localization property of a random chain changing under the influence of a constant electric 
field has been studied. We have adopted the multifractal scaling formalism to explore the possible 
localization behavior in the system. We observe that the possible localization behavior with the 
increase of electric field is not systematic and shows strong instabilities associated with the local 
probability variation over the length of the chain. The multifractal scaling study captures the local- 
ization aspects along with strong instability when the electric field is changed by infinitesimal steps 
for a reasonably large system size. 
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Electronic states are exponentially localized in one di- 
mensional (hereafter ID) random chains and the enve- 
lope of the wave function, (j){x) ~ exp {—ax) for x — > oo 
where a is the inverse of localization length . This 
localization nature of electronic states can be changed 
through application of a constant electric field and as a 
result electronic states exhibit some kind of different lo- 
calized nature over the sample size. The problem of elec- 
tronic states and localization in a random chain in the 
presence of a constant electric field is still a matter of con- 
troversy and is yet to be fully understood In the 
past, the possibility of existence of non-exponential local- 
ization or localization/delocalization transition in electri- 
fied chain has been addressed HJ^. The deviation from 
exponential localization has been also claimed through 
the numerical study of electronic transmittance ^ . 

However, the localization mechanism can be under- 
stood more rigorously within the multifractal scaling 
analysis of the electronic wave functions without any a 
priori assumption of exponential localization nature and 
the existence of localization length pO| , p^ -p^ . The mul- 
tifractal scaling formalism has been invoked in the 
recent past to analyze the nature of electronic states in 
the vicinity of the mobility edge [p"2|jl3| , and also for 
characterization of critical nature of electronic states in 
ID Fibonacci quasiperiodic systems ]l5| , [T6t . Generally 
speaking, in all of the above examples, wave functions 
exhibit a rather involved oscillatory behavior displaying 
strong fluctuations. As a consequence the notion of en- 
velope wave function or lyapunov exponent which has 
been successful in studying both the extended and ex- 
ponentially localized states, is no longer suitable to the 
states in examples above. On the other hand multifrac- 
tal scaling formalism has been found to be very useful to 
characterize the spatial fluctuating pattern of the wave 



functions which are neither Bloch like homogeneously ex- 
tended state nor the exponentially decaying one. The 
same scaling analysis has been successful to extrapolate 
in the extreme limits also. So, it naturally appears that 
the possible delocalization behavior that we are going to 
address in this Letter can also be understood through the 
same scaling analysis. 

The aim of this Letter is to report our investigation 
on how does localization nature is influenced by switch- 
ing on a constant electric field through numerical study 
of electronic wave functions for a reasonably large array 
of (5-function random potentials having a bi-modal dis- 
tribution. The choice of this type of potential has been 
justified in the past both from its experimental relevances 
as well as from the pure academic interest. We start with 
the Schrodinger's equation for the electrons in a random 
chain in presence of a constant electric field : 
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^{x) = E^{x) (1) 



where the units arc such that [Ti = 2me = 1) mg is the 
effective mass of electron. The electric field induced force 
F is expressed in unit of {^ 2m ) ' ^ lattice spac- 

ing and Vn is the strength of the n-th potential barrier 
(taking the value Va or Vb randomly), F is the prod- 
uct of electric field by the electronic charge. The lattice 
constant a is taken unity throughout this calculation. 

One can map the above Eq.(|l|) to a finite difference 
equation by approximating the potential Fx by a step 
function in-between the i5-functions. Within this approx- 
imation, the solutions in-between the (5-function poten- 
tials are now plane waves instead of Airy functions. The 
corresponding Poincare map P,pUr^ is: 
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= A„*„ + B„*„_l 

The coefficients An and i?„ are given by § 
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with fc„ = (£; + nF)i/2 and v|/„ = = n). Now 

in order to solve the equation iteratively for a reason- 
ably large system size one can consider the initial values 
for = eyi-p{-iE^/^a) and *2 = {-2iE''-/^a), E being 
the incident electron energy before it reaches the region 
where the electric field is applied. 

The transmittance corresponding to the array of ran- 
dom (5-function potential is given by B 



exp(2ifci) — 1 
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with 



k = E^/^ and ki = (E + FLf''^ 



and L = Na. 

We now analyze the pattern of local probability density 
|\1/„P along the chain through its multifractal scaling re- 
lation of a and /(a), where a stands for the scaling expo- 
nent and f{ot), the corresponding distribution function. 
We have used the mathematical prescription suggested 
by Chabbra and Jensen ||l^ for its simplicity and success 
in correctly evaluating the quantities a and f{a) directly 
through the normalized measure without any numerical 
instability. Let us define the required normalized mea- 
sure in our study by 
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where the scaling behavior of Pi ~ iV~"' for N ^ oo. 

According to Chabbra and Jensen if we define the g-th 
moment of the probability measure Pi by fj,i (g, N) where 



P? 



then a complete characterization of the fractal singular- 
ities can be made in terms of ^i(q^N). The expression 
for the distribution function of scaling exponent a can 
be written as: 



1 



N 



f(a) = hm — , / 
° i=i 



fi^iq,N) log fi,{q,N) (5) 



and the corresponding singularity strength of the mea- 
sure is obtained by 
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One can infer on the nature of electronic states for 
large N based on the following observation: 



1. Extended nature : an 
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2. Localized nature : amin 
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This property is usually manifested by rectangular 
two-hump form of /(a) curve with a sparse distri- 
bution of points in between. 

3. Critical nature : a vs f{a) curves closely overlap 
on one another with the increase of system size. 

4. Power-law nature : The right portion of a — /(a) 
curve deviates slowly from one another with the 
increase of system size in contrast to the strong 
deviation as seen in the exponential decay. 

We now define a simple way the degree of localization. 
We consider Aq, = {amax — cimin) as a measure of degree 
of localization. This can distinguish clearly an extended 
state from a localized state for a sufficiently large system 
size N. Also, one can investigate the change in the na- 
ture of states brought about due to the change of some 
external parameters, e. g., electric field. 

In figure (1) we have shown the spatial pattern of local 
probability variation along the chain for a localized state 
in both zero and finite electric field. The upper curve ex- 
hibits the delocalized pattern in the presence of electric 
field F — 1.25 X 10~^ unit. The corresponding localized 
and delocalized behavior of the electronic transmittance 
data have been been presented for both the zero as well 
as for a finite electric field in figure (2). In figure (3), 
we have shown (a — f{a)) plots for both the zero and 
the finite field value 2 x 10^^. In the zero field case the 
plot shows in f{a) a two-hump form corresponding to 
the exponential localization. This is due to the fact that 
(a — f{a)) spectrum is densely populated on the extreme 
left and on the right but the region in-between is very 
sparse. On the other hand, for finite electric field, the 
(a — f{a)) data gather in a relatively narrow region on 
the curve having a convex shape. This indicates that 
the change in localization behavior shows an overall de- 
localization trend which is exhibited through the more 
spatial extension of the state over the sample and hence 
resulting in the contraction of the a — f{a) spectrum. 

Next we investigate further whether this delocalization 
pattern can change systematically as we change the elec- 
tric field through small steps of the order of 10~* unit for 
a sufficiently large system size. 

In figure (4) we have shown the plot of degree of lo- 
calization Aq, with the electric field for two large system 
sizes, N = 10^ (upper curve) and = 2 x 10^ (lower 
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curve). In both the plots we sec that Aq, exhibits strong 
instabihty throughout the whole regime of electric field 
from a very low to as high as 10~^ unit. Wc have ob- 
served that for both the two system sizes, the order of 
deviation in from its previous step is of the order of 
unity whereas the value of A^ itself is of the same order. 
Also the fluctuating pattern of A„ with F for the both 
the two system sizes is almost the same. In figure (5) we 
have shown again the variation of the degree of localiza- 
tion, i. e., Aa with the electric field for a different set 
(cf. Figure 4) of potential parameters. Here the order of 
fluctuations in Aq, corresponding to two large systems, 
i. e., 75 X 10^ and 1.7 x 10^ number of atoms have been 
presented. The order of fluctuations appear to be quite 
significant and it is nearly unity in the both the cases 
over a wide region of the electric field as shown in the 
figure through covering the fluctuating zone between the 
horizontal lines. 

We think this kind of instability has its intrinsic origin 
in the restructuring of the different states in a compli- 
cated manner and each of them is highly sensitive even 
for an infinitesimal change in the applied electric field. 
This can be also understood as due to the competing na- 
ture of the potential Fn and the disordered (5-potentials 
in the large system size. However, if one neglects this 
fluctuation through some brute force methods, it shows 
only an apparent simple delocalization effects due to the 
increase of electric field. 

In conclusion, we have shown that the change of local- 
ization aspects due to the increase of a uniform electric 
field is not simply an overall delocalization but the local- 
ization property is very sensitive with respect to an in- 
finitesimal change of electric field giving rise to a strong 
instability in the degree of localization. This instability 
aspect in the localization/delocalization is present for all 
reasonably large lengths for an appropriate set of param- 
eters and is due to the combined effects of disorder po- 
tential and the electric field induced linear potential. At 
sufficiently large length scale the states change depending 
upon the restructuring of the spectrum from its previous 
form and hence the localization property of a particular 
state in a given field will change drastically giving rise 
to a state of modified localized nature. P. Biswas would 
like to thank the Council of Scientific and Industrial Re- 
search (CSIR) for financial assistance in the form of a 
senior research fellowship. 
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FIG. 1. Local probability (in log-scale) vs system size 

{N) in unit of 10* number of sites. The parameters arc : 
Va = 1.25 units, Vb = 1.4 units, c = 0.01, AT = 4 x lO'^ for 
electric field F = 2.0 x 10"^ (upper) unit and iV = 5.6 x 10* 
for zero field (lower) with the energy of the incoming electron 
E = 1.4890 units in both cases. The data in the upper curve 
has been given a constant shift of 20 units in the log-scale for 
convenience of comparison. 
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FIG. 2. Transmittance T (in log-scale) vs size [N) for the 
set of parameters as in figure 1. The plot is made up to a 
length of 4 x 10® number of sites for F = 2.0 x 10~® (upper 
curve) and up to 5 x 10* number of sites in the zero field case 
(lower curve). 
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FIG. 3. a vs /(a) plot for finite {F ■■ 
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: 2 X 10"® unit) as well 
as zero electric field for the local probability variation 
over the chain of length iV = 4 x 10^ . The other parameters 
are the same as in figure (1). 
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FIG. 5. The degree of localization (Ac) vs electric field 

(F) for a different set of potential parameters. The thick 
and the dashed line corresponds to system size 7.5 X lO"* and 
1.7 X 10^ with other parameters, Va = 5.15 units, Vb = 5.34 
units, and c = 0.05, and energy E = 5.6025 units. 
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FlinunitoflO'l-e)) 

FIG. 4. The degree of localization (Aq) vs electric field 
(F) for the two different system sizes 10® (thick line) and 
2 X 10® (dashed line). The parameters are Va = 1-25 units, 
Vb = 1.4 units, c = 0.01, and for incoming energy E = 1.4890 
units. Once again, for comparison, we plot the negative of 
for system size 2 x 10® (dashed line). 
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